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ABSTRACT 

The debris dispersion model (which dispersion model?) 
is so heterogeneous and interrelated with various 
factors, 3D graphics combined with physical models are 
useful in understanding the complexity of launch and 
range operations. Modeling and simulation in this area 
mainly focuses on orbital dynamics and range safety 
concepts, Including destruct limits, telemetry and 
tracking, and population risk. Particle explosion modeling 
is the process of simulating an explosion by breaking the 
rocket into many pieces. The particles are scattered 
throughout their motion using the laws of physics 
eventually coming to rest. The size of the foot print 
explains the type of explosion and distribution of the 
particles. The shuttle launch and range operations in 
this paper are discussed based on the operations of the 
Kennedy Space Center, Florida, USA. Java 3D graphics 
provides geometric and visual content with suitable 
modeling behaviors of Shuttle launches. 

INTRODUCTION 

Over the past several years there has been a steady 
increase in space launch operations. NASA is 
developing operational concepts that will seamlessly 
integrate air and space launch/reentry operations while 
ensuring that the highest levels of safety, security and 
efficiency are maintained. For launches at Kennedy 
Space center and Cape Canaveral, range safety is 
maintained from the range operations control center 
iocated at Cape Canaveral Air force station 
[1,12,13,14,15].. Debris dispersion is a key safety 
consideration in the launch decision process, and 
understanding the vehicle’s design characteristics, 
atmospheric conditions, and other factors that can affect 
a hazard impact area is a critical issue for range 
operators. To evaluate and mitigate risks within the 
range, range operators use the Common real time debris 
foot print system. In general, a debris dispersion model 
computes the velocity of fragments, size distribution, and 
the impact probabilities and risk associated with debris 
resulting from a catastrophic failure or intentional 


destruction of a vehicle [7, 1 7]. Presently the eastern 
range uses the CRTF (Common Real-Time Debris 
Footprint) [5] program operating inside the RRAT (Range 
Risk Analysis Tool) to simulate debris dispersion and 
perform the risk analysis. In the next section, the basic 
physics of falling debris are discussed. 

MECHANICS OF DEBRIS FALL AND 
DISPERSION 

A piece of debris has an initial state vector that is defined 
by a position and velocity vector. The initial state vector 
can be perturbed from by an explosion that imparts a 
velocity and a consequential adjustment to the velocity 
vector. There is no adjustment to the initial position 
because the velocity is added impulsively.. The 
dominant parameter in the trajectory computation is the 
ballistic coefficient [10]; a measure of an object’s ability 
to overcome air resistance. Objects with low weight to 
drag ratio (low ballistic coefficient) fall much more slowly 
than objects with a high weight to drag ratio (high ballistic 
coefficient). Hence, fragments with low ballistic 
coefficients will tend to travel further given an initial 
horizontal velocity. In this situation wnd becomes a major 
factor affecting the fall of the debris. 
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Figure 1 Generic Debris Dispersion Mode! 


The dispersion of debris impact location associated with 
a given initial state vector is illustrated in figure 1 . The 
primary sources of debris dispersion are ballistic 
coefficient, wind and velocity perturbation. Wind data is 



obtained from Rawinsonde data [8] with upper air 
* . characteristics. The following equations provide 
dispersion of particles w.r.t x,y,z position of translation in 
3-D space. 

X= X p +V x *At+W x *At+ V x *C x ‘At 

Y= Y p +V/At+W y *At+ V y *C y *At 

Z= Z p +V z *At+W z *At+ V z *C z *At-0.5*9.81 *At*At 

Where p indicates the position where explosion occurs 
(you state above that the explosion is an “impulsive” 
process, therefore it doesn't have a “start” or “end”), W 
indicates the wind local velocity and direction impact on 
the individual particles (Question: does this include the 
ballistic coefficient)?, C expresses the coefficient of 
explosion and explostion type, which is directly related to 
the dispersion effect. At is the time increment of the 
simulation.. The simulation treats each piece of debris 
as a unique projectile and uses the above projectile 
equations to calculate Its state at each time step. This 


simulation runs until all the particles hit the ground. 
(Question: do you consider spherical geometric effects?) 
The air resistance and elasticity effecis are also added to 
particles for a reality effect.. In the model, the launch 
pad at Kennedy Space Center is considered as point of 
origin (0,0,0). The reference co-ordinate system X refers 
to east, Y pointing to the true north and Z is the altitude. 
All measured units in the model are presented in SI 
units. With the above basic physical theory and 
concepts, a JAVA 3D visualization for debris dispersion 
is developed and explained in following sections. 

JAVA 3D MODEL 

This application has been developed as an Applet which 
incorporates physical model in Java which computes 
rocket characteristics like velocity, altitude, downrange 
etc. and Java 3D dispersion model [16]. The graphical 
user interface and output of dispersion model is shown in 
figure2. 



Figure 2. Graphical user interface of debris dispersion model 




The left hand side of the the applet accepts launch and 
rocket characteristics that allow computation of rocket 
velocity and down range distance. The results are 
displayed as a 3-D animation showing the rocket’s 
trajectory and destruct lines [2,4], The user is allowed to 
specify up to a three stage rocket There are three stages 
allowed, but the visualization effect of debris dispersion 
model applies for the first two stages of launch. 

MULTISTAGE TRAJECTORY ANALYSIS 

The general characteristics of launch details like latitude 
of launch pad (degrees), starting altitude (m), starting 
velocity (m/s) are provided in the input form in figure 2. 
The specific stage details like empty mass ratio for a 
given stage, exhaust velocity (m/s), liftoff acceleration 
(m/s 2 ), payload ratio and initial mass (kg) are required to 
compute trajectory details. Airspeed is the speed of the 
spacecraft relative to the launch site. Airspeed along 
with the deceleration; which is derived from the density 
and frontal area of the spacecraft, is used to calculate 
the air resistance. 

air resistance = deceleration * airspeed 2 

When a spacecraft reaches low earth orbit, it is generally 
at an altitude of one hundred and forty kilometers A high 
starting altitude is desirable because the air pressure 
and density decrease as the altitude increases. A lower 
air density means the aerodynamic drag will be lower 
and a lower air pressure means the exhaust velocity will 
be higher. Altitude is used to calculate air density and 
orbital velocity. 

air density = sea level density * exp( - altitude / column height ) 
orbital = ground orbital * sqrt( radius / ( radius + altitude ) ) 

Column height is defined as the average column height 
of the atmosphere; the height rise after which the 
pressure falls by a factor of e. This is calculated using 
gravity, molecular weight and temperature. This is then 
used to derive the air density as the spacecraft rises. 
The Radius refers to earth’s radius. 

Column height = 8314 ( J * K / kmol ) * temperature / gravity / 
molecular weight 

Air density = sea level density * exp( - altitude / column height ) 

Downrange is the horizontal distance traveled by the 
spacecraft. When a spacecraft reaches low earth orbit, it 
is generally downrange about nine hundred kilometers , 
however it is dependant on several launch conditions. 
The downrange when reaching orbit will decrease as the 
starting velocity increases. The downrange when 
reaching orbit will increase as the acceleration 
decreases. For example, a low acceleration space plane 
will be far downrange when reaching orbit, while a high 
thrust rocket will be relatively close to the launch position. 


The velocity and downrange distance are computed 
using above equations, which form the initial velocity, 
altitude and down range distance from launch pad when 
explosion occurs. In this model, the number of particles 
is assumed as 500 random size pieces. The user can 
alter any assumptions of payload, vehicle types and 
launch pad scenarios in Kennedy Space Center. 

VISUALIZATION 

The scene graph consists of superstructure components, 
a VirtualUniverse object and a Locale object, and a set of 
branch graphs which is displayed via 3D canvas [4,6]. 
The simulation is done via Java 3D behaviors. There are 
three behaviors which correspond to orbit, explosion and 
dispersion of particles. Behaviors allow specifying an 
action based on a particular set of events called the 
WakeupCondition. A WakeupCondition is a combination 
of another Java 3D object called a WakeupCriterion, 
which covers a wide range of activities like keyboard 
stroke, frames, waiting, etc. The spaceHit behavior 
enacts a key board strike, which removes the shape of 
the rocket from the scene graph, and adds a 
‘createexplosiongrapH branch graph which displays a 
representation of an explosion. The next key stroke 
(what was the first key stroke?) invokes the explosion 
behavior, which is a specialized interpolator. This 
interpolator describes the dispersion of the particles from 
the point of explosion and scattering around the path of 
dispersion to the ground continuously. The transition 
from the starting to ending value is interpolated 
continuously from the value 0-1.0. The Alpha class is a 
simple object that help Java 3D worlds come alive. 
Alpha object generates a single ramp between 0 and 1 .0. 
There is a random number generator allocating which 
trajectory rocket should be launched. If rocket entered 
into an abnormal trajectory or crossed destruct limit lines, 
spaceHit behavior is invoked to kill the rocket. The 
VirtualUniverse for debris dispersion model is shown in 
figure 3. 

If the rocket orbits in a perfect trajectory and 
corresponding Alpha reaches 1.0, the satellite orbit is 
activated in another universe (??). The trajectories are 
constructed by Bezier curves [11], The origin, end point 
and another two points which represent between the 
start and end points are called control points. The 
curves are constructed as a sequence of cubic 
segments, rather than linear ones. The entire curve is 
contained in the quadrilateral whose comers are the four 
given points (their convex hull). Constructing these 
curves is very efficient due to the employment of a 
simple recursive algorithm that reduces the computation 
of points to a simple “division-by-two” operation. Seven 
trajectories are constructed and dotted lines represent 
destruct limit lines. The bottom texture of the box has 
the appearance of Kennedy Space Center with launch 
pad locations in figure 4. 
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Figure 3 Scene graph for Shuttle launch and debris dispersion model 


Figure 4 Debris dispersion model visualization 












The debris dispersion is displayed in figure 4 showing 
that all debris particles fall into tie elliptical impact zone 
(orange cluster of particles in right hand side over Florida 
and Atlantic Ocean). The debris dispersion model 
obtains Rawinsonde data [8] from an online weather 
server and incorporates wind effects into debris 
dispersion model computations shown above. The three 
dimensional visualization provides the debris scattering 
impact zone the potentially affected areas. Based on 
debris impact, human health risk assessment and 
expected casualty can be derived [3,9,19]. 


CONCLUSION 

Earlier debris dispersion model developed are based on 
two-level full factorial analysis on the CRTF [5] debris 
dispersion using experimental design techniques [18]. 
From the above model, the results show that three 
factors have statistical significance; ballistic coefficient of 
debris fragments (C), failure altitude (F) and wind speed 
and direction (W). In the model presented, the Java 3D 
model shows a simple orbital and debris dispersion 
model based on physics with necessary parameters and 
their combined interactions in each fragments. The 
model will be used by Range Safety Officers as a 
simulation tool during decision making process.. Future 
research will focus on an updated dispersion model 
combined with all available parameters using Java 3D. 
Future model make a decision to destroy Shuttle when it 
violates flight launch rules with built in intelligent systems. 
Developing the application in Java 3D aids in its 
distribution and usefulness across all NASA centers. 

The validation of models is based on telemetry and radar 
data to monitor orbital dynamics and any debris 
dispersion. Outside the Shuttle program, it can be used 
for missile tracking or any trajectory analysis. The 
visualizations are self-explanatory and can be interpreted 
easily by user. This visualization tool can be analyzed 
collaboratively using web. The applications of this tool 
can be used for studying any catastrophe in launch and 
range operations in lunar and mars surface. This model 
can also be integrated into onboard systems for 
telemetry tracking and onboard destruct mechanisms 
based on flight rules. 
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